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The processes 77 — ujcj), (fxf), and luuj are measured using an 870 fb~^ data sample collected with 
the Belle detector at the KEKB asymmetric-energy e"'"e~ collider. Production of vector meson pairs 
is clearly observed and their cross sections are measured for masses that range from threshold to 
4.0 GeV. In addition to signals from well established spin zero and spin two charmonium states, 
there are clear resonant structures below charmonium threshold, which have not been previously 
observed. We report a spin-parity analysis for the new structures and determine the products of 
the rjc, XcO, and Xc2 two-photon decay widths and branching fractions. 



PACS numbers: 14.40.-n, 13.25.Gv, 13.25.Jx, 13.66.Bc 

A plethora of states, especially many new charmo- 
nium or charmonium-like states (the so called "XYZ 
particles"), that are not easily accommodated vifithin the 
quark model picture of hadrons have been observed [ij. 
Recently a clear signal for a new state X{3915) — )• 
ujJ/ip i2] and evidence for another state X(4350) 
(jjJ/ip [3] have been reported, thereby introducing new 
puzzles to charmonium or charmonium-like spectroscopy. 



Since these states couple to a J/^p and a light mass vector, 
some authors have suggested that they are good candi- 
dates for molecular or tetraquark states [l|. 

It is natural to extend the above theoretical picture 
to similar states coupling to uj(j), since the only difference 
between such states and the X(3915) ^ or X(4350) fS\ is 
the replacement of the cc pair with a pair of light quarks. 
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A good neutral cluster is reconstructed as a photon if 
its electromagnetic calorimeter (ECL) shower does not 
match the extrapolation of any charged track and its 
energy is greater than 50 MeV. The tt" candidates are 
reconstructed from pairs of photons with invariant mass 
within 15 MeV/c^ of the tt" nominal mass. Here the 
7r° mass resolution is about 6 MeV/c^ from MC simu- 
lation. A mass-constrained kinematic fit is applied to 
the selected 7r° candidate and < 10 is required. For 
77 — ^ WW, the energies of the photons from 7r° decays 
are further required to be greater than 75 MeV in the 
the pOp° threshold that is absent in 77 p+ p- [|. The endcap ECL region (cos 61^ < -0.65) to suppress back 



States coupling to ww or although not as exotic as 
those that decay into W(^, which have two pairs of light 
quarks with different flavor, could also provide informa- 
tion on the classification of the low-lying states coupled 
to pairs of light vector mesons. 

Experimental studies of 77 VV (V — p, w, 0, K*) 
began in 1980 with the measurement of 77 p'^p^ [3], 
and later 77 — > p'^ p~ 0. A number of theoretical mod- 
els, such as q'^q^ tetraquark statesjS], Regge exchange j7[], 
and an s-channel p^p" resonance [8|, were proposed to ex- 
plain the large cross section observed in 77 p^p^ near 



77 — >• wf/) and WW pro cesses were studied by the ARGUS 
Collaboration [13, liil with very limited statistics, while 
77 — >■ 00 has never been measured below the charmo- 
nium mass region. 

In this Letter, we report measurements of the cross 
sections for 77 — VV, where VV = uj(f), (fxj) and ww, 
as well as observations of new resonant structures be- 
low charmonium threshold. The results are based on an 
analysis of an 870 fb~^ data sample taken at or near 
the T(nS) {n — 1, ...,5) resonances with the Belle detec- 
tor 'l2| operating at the KEKB asymmetric-energy e~^e~ 
cohider [ll. 

The Belle detector is described in detail elsewhere 



We use the program treps ;14j to generate signal Monte 
Carlo (MC) events and determine experimental efficien- 
cies and luminosities. In this generator, the two-photon 
luminosity function is calculated and events are gen- 
erated at a specified fixed two-photon center-of-mass 
energy (W^^) using the equivalent photon approxima- 
tion [15| . The efficiencies for the detection of 77 —5- X 



W(/ 

and 77 — > X WW — > 2(7r+7r~7r°) are determined as- 
suming different spin-parity ( J^) assignments and a zero 
intrinsic width for the X. 

We require four reconstructed charged tracks with zero 
net charge. For these tracks, the impact parameters 
perpendicular to and along the beam direction with re- 
spect to the interaction point are required to be less 
than 0.5 cm and 4 cm, respectively, and the transverse 
momentum in the laboratory frame is restricted to be 
higher than 0.1 GeV/c. For each charged track, infor- 
mation from different detector subsystems is combined 
to form a likelihood Ci for each particle species [l^. A 



track with TZk 



> 0.6 is identified as a kaon. 



while a track with TZk < 0.4 is treated as a pion. With 
this selection, the kaon (pion) identification efficiency is 
about 97% (98%), while 0.4% (1.0%) of kaons (pious) are 
misidentified as pious (kaons). A similar likelihood ratio 
is formed for electron identification [l3|- Photon con- 
version backgrounds are removed if any charged track in 
an event is identified as electron or positron {TZe > 0.9). 
For 77 (j)(f), we require that only three of the charged 
tracks be identified as kaons. 



ground with misreconstructed photons. When there are 
more than two tt'^ candidates in an event, the pair with 
the smallest sum from the mass constraint is retained. 
To suppress backgrounds with extra neutral clusters in 
the w0 and ww modes, events are removed if there are 
additional photons with energy greater than 160 MeV. 

We define the w signal region as 0.762 GeV/c^ < 
M(7r+7r-7r°) < 0.802 GeV/c^, and the w mass sidebands 
region as 0.702 GeV/c^ < M(7r+7r-7rO) < 0.742 GeV/c^ 
or 0.822 GeV/c2 < M(7r+7r-7rO) < 0.862 GeV/c^, which 
is twice as wide as the signal region. The signal 
region is defined as 1.012 GeV/c^ < M{K+K-) < 
1.027 GeV/c^, and its sideband regions are defined 
as 0.99 GeV/c2 < M{K+K-) < 1.005 GeV/c^ or 
1.034 GeV/c^ < M{K+K-) < 1.049 GeV/c^. The 
(j) sidebands are also twice as wide as the signal re- 
gion. For the two possible combinations of (j)(f> in the 
2{K~^K~) final state, the one with the smallest 6min — 
^{M{K+K-)i ~ m^Y + {M{K+K-)2 - m^Y is cho- 
sen. For the four possible combinations of ww, only one 
combination from a true signal can survive after event 
selection. 

The magnitude of the vector sum of the final particles' 
transverse momenta in the e^e~ center-of-mass (CM.) 
frame, |X]-Pt*|j which approximates the transverse mo- 
mentum of the two-photon-collision system, is used as 
a discriminating variable to separate signal from back- 
ground. The signal tends to accumulate at small | X] -F**! 
values while the non-77 background is distributed over 
a wider range. We obtain the number of VV events in 
each VV invariant mass bin by fitting the | X] ^t* I distri- 
bution between zero and 0.9 GeV/c. The signal shape 
is from MC simulation of the signal mode and the back- 
ground shape is parameterized as a second-order Cheby- 
chev polynomial. In order to control the background 
shape, we restrict the coefficients of the background poly- 
nomials in nearby invariant mass bins to vary smoothly 
along parabolas. The parameters of these parabolas are 
determined from fits to the coefficients obtained from fits 
to the I X] ^t* I distribution in each VV invariant mass bin. 
The resulting VV invariant mass distributions are shown 
in Fig. m 

In Fig. [1] there are some obvious structures in the low 
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FIG. 1: The (a) a;0, (b) (p(j) and (c) ljuj invariant mass dis- 
tributions obtained by fitting tlie | -Pt* | distribution in each 
VV mass bin. The shaded histograms are from the corre- 
sponding normalized sidebands, which will be subtracted in 
calculating the final cross sections. 



VV invariant mass region. Two-dimensional (2D) an- 
gular distributions are investigated to obtain the of 
the structures. In the process 77 VV, five angles are 
kinematically independent. Among the possible variable 
sets, we choose z, z*, z**, cj)* , and (j)** . Using ojcj) as an 
example, z is the cosine of the scattering polar angle of 
(j) in the 77 CM. system; z* and (f>* are the cosine of the 
helicity angle of K'^ in the (j) decays and the azimuthal 
angle defined in the rest frame with respect to the 
'Y^ LUip scattering plane; z** and </>** are the cosine of 
the helicity angle of normal direction to the decay plane 
of the uj — Tr+TT^Tr*^ and the azimuthal angle defined in 



the UJ rest frame. We use the transversity angle (0t) and 
polar- angle product (He) variables to analyze the angular 
distributions. They are defined as (px — \(t)* + 0**|/27r, 

iig = [i-{z*ni'-iz**n 

We obtain the number of signal events by fitting the 
^ Pj* I distribution in each pT and Ilg bin in the 2D 
space , which is divided into 4 x 4 , 5x5, and 1 x 1 bins for 
4>4>, and ujuj, respectively, for M{VV) < 2.8 GeV/c^. 
The obtained data are fitted with the signal shapes from 
MC-simulated samples with different assumptions 
(0+,0-,2+,2-). The fitted results are: (1) for top: a 
mixture of 0+ (S'-wave) and 2+ (S'-wave) describes data 
with /ndf — 0.9 {ndf is the number of degrees of free- 
dom); (2) for (fxp: a mixture of 0+ (S'-wave) and 2~ {P- 
wave) describes data with /ndf — 1.3; and (3) for ujlu: 
a mixture of 0+ (S'-wave) and 2+ (S-wave) describes data 
with 'X^ jndj = 1.3. The contributions from other are 
found to be small and thus neglected. 



The cross section cr^^_j.yy (W-y^) is calculated from 



Ar 



dh 



dW. 



(1) 



' 77 



1 dTj ^ 

where — — 



dW^ 



is the differential luminosity of the two- 
photon collision, and e is the efficiency. Here ISW^^ is 
the bin width and An is the number of events in the 
AW-y^ bin. 

The 77 VV cross sections are shown in Fig. [2] For 
the processes 77 ucj) and (jxj), the cross sections are 
measured in the CM. angular range jcosS*! < 0.8 since 
there are no detected events beyond this limit. In cal- 
culating the cross sections, J^- weighted efficiency curves 
are used based on the fitted angular distribution results. 
The effect of the mass resolution (which is a few MeV/c^) 
is neglected since it is much smaller than the bin width. 
The fraction of cross sections for different values as 
a function of M{yV) is also shown in Fig. [21 We con- 
clude that there are at least two different components 
( J = and J = 2) in each of the three final states. 

The inset also shows the distribution of the cross sec- 
tion on a semi-logarithmic scale, where, in the high en- 
ergy region, we fit the dependence of the cross sec- 
tion. The solid curves are the fitted results; the fit gives 
n = 7.2 ± 0.6, 8.4 ± 1.1, and 9.1 ± 0.6 for the ww, 
and (jxj) modes, respectively. These results are consis- 
tent with the predictions from pQCD and handbag mod- 
els [l8|. and similar to previous measurements in other 
modes [l9|. 

There are several sources of systematic error for the 
cross section measurements. The particle identification 
uncertainties are 1.5% for each kaon, 1.2% for each pion. 
An uncertainty in the tracking efficiency for tracks with 
transverse momenta greater than 200 MeV/c is about 
0.35% per track. For tracks with low transverse mo- 
menta, the uncertainty is estimated by studying the 
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FIG. 2: The cross sections of 77 — >■ u}(j> (a), (jxj) (b), and ljuj 
(c) are shown as points with error bars. The fraction contri- 
butions for different values as a function of M{VV) are 
shown as the points and squares with error bars. For the 
processes 77 — >■ ujcj) and (fxp, the cross sections are measured 
in the CM. angular range |cosS*| < 0.8. The error bars are 
statistical only; there are overall systematic errors of 15%, 
11% and 13% for ujcp, cj>(j> and uui, respectively. The inset 
also shows the cross section on a semi-logarithmic scale. In 
the high energy region, the solid curve shows a fit to a W^^ 
dependence for the cross section after the significant charmo- 
nium contributions {rfc, XcO and Xc2) were excluded. 
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control samples. A 



momentum-weighted systematic error is taken for each 
track. The efficiency uncertainties associated with the lo 
and (j) mass requirements are almost independent of the 
VV mass, and are estimated to be 1.9% and 1.6%, respec- 
tively. The statistical error in the MC samples is about 



0.5%. The accuracy of the two-photon luminosity func- 
tion calculated with the treps generator is estimated 
to be about 5% including the error from neglecting ra- 
diative corrections (2%), the uncertainty from the form 
factor effect (2%), and the uncertainty in the total inte- 
grated luminosity (1.4%) IJ]. The trigger efficiency for 
four charged track events is rather high because of the 
redundancy of the Belle ffist level multi-track trigger. 
According to the MC simulation, the trigger efficiencies 
for signal increase with increasing VV mass. The uncer- 
tainty of the trigger simulation is smaller than 5% [20| . 
The preselection efficiency for the final states has little 
dependence on the VV invariant mass, with an uncer- 
tainty that is smaller than 1% for uj(f), 4% for (fx}) and 
2.5% for LOU. From Ref. [2l[, the uncertainty in the world 
average values for B{(t> K'^K~) is 1.1% and that for 
B(ijj 7r"'"7r~7r°) is 0.8%. The uncertainty in the fitted 
yield for the signal is estimated by varying the order of 
the background polynomial and fit range. The uncer- 
tainty on the I X) -P** I resolution is estimated by changing 
the MC signal resolution by ±10%. The uncertainty on 
the weighted efficiency curve is estimated by changing 
the fitted ratio of the components by ztltr. Assuming 
that all of these systematic error sources are independent, 
we add them in quadrature to obtain the total systematic 
errors, which are 15%, 11% and 13% for cjt/), (fxf) and ujuj, 
respectively. 

For VV invariant masses above 2.8 GeV/c^, we mea- 
sure the production rate of charmonium states. For 
these measurements, | X) -^t* I required to be less than 
0.1 GeV/c in order to reduce backgrounds from non- two- 
photon-processes and two-photon-processes with extra 
particles other than (/) or a; in the final state. 

Figure [3] shows the VV invariant mass distributions. 
Clear ric, Xco and Xc2 4>4'i and ric signals are ev- 

ident. The VV mass distributions are fitted with three 
incoherent Breit-Wigner functions convoluted with a cor- 
responding double Gaussian resolution function as the r?c, 
Xco and Xc2 signal shapes, and a second-order Chebychev 
polynomial as the background shape. The shape of the 
double Gaussian resolution function is obtained from MC 
simulation. From the fits, we obtain yields of 386 ± 31, 
56±11 and 89±11 events in and 85±29, 19±10 

and 16 ± 7 events in ww for rjc, XcO and Xc2 signal events, 
respectively. The charmonium yields are consistent with 
zero in the ujcj) mode. Bayesian upper limits on signal 
yields are estimated to be 7.9, 4.3 and 2.4 for rjc, XcO 
and Xc2 i^4>j and 35 and 28 for XcO and Xc2 
respectively, at the 90% confidence level (2l| . 

The product of the two-photon decay width and 
branching fraction T^^B{X VV) (or the upper lim- 
its in case the signal is insignificant) for rjc, XcO and Xc2 
are listed in Table IH A similar systematic error estima- 
tion is performed together with the uncertainties in the 
resonance parameters results in the total systematic er- 
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rors of 13%, 11%, and 11% for r^^{R)B{R uj^); 7.9%, 
8.0%, and 7.2% for r^^{R)B{R -> 00); and 11%, 10%, 
and 9.1% for T^^{R)B{R — > ojlu), for R ~ rjc, Xco and Xc2i 
respectively. For the upper limit determinations, the ef- 
ficiencies have been lowered by a factor of 1 — Csys in 
order to obtain conservative values. The measurements 
o{Tj^B{X (jxp) for r]c, Yr n an d Xc2 are consistent with 
previously published results [20| with improved precision. 
The value of T~f^B{X 4'4>) for r/c, Yrn a,nd Xc2 obtained 
in this work supersedes that in Ref. [20||. All the other 
results are first measurements. 
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FIG. 3: The invariant mass distributions of (a) u)(f), (b) 4>ij>, 
and (c) uu) combinations in the charmonium mass region. The 
points with error bars are data, and the solid curves are the 
best fits. 

In summary, we present a search for exotic states in 
two-photon processes 77 uj(f), (fxp and louj. Clear 



(jjcf), (j)(j>, and ujoj production is observed, and cross sec- 
tions are measured. We observe clear structures at 
M(a;0) - 2.2 GeV/c^, M(00) - 2.35 GeV/c^, and 
M{ujuj) ^ 2.0 GeV/c^. While there are substantial spin- 
zero components in all three modes, there are also spin- 
two components near threshold. The measured cross sec- 
tions for 77 — > cjt/) and 00 are comparable in the mass re- 
gion between 2.0 GeV/c^ and 2.5 GeV/c^, which disagree 
with the intuitive expectation that the former is harder to 
produce than the latter since the quark flavors are differ- 
ent in w and 0. The cross sections for 77 uicj) are much 
lower than the prediction of the q'^q'^ tetraquark model Q 
of 1 nb, while the resonant structure in the 77 (fxj) 
mode is nearly at the predicted position. However, the 
00 cross-section is an order of magnitude lower than the 
expectation in the tetraquark model. On the other hand, 
the t-channel factorization model [23j] predicted that the 
00 cross sections vary between 0.001 nb and 0.05 nb in 
the mass region of 2.0 GeV/c^ to 5.0 GeV/c^, which are 
much lower than the experimental data. For 77 — >■ ww, 
the t-channel factorization model Q predicted a broad 
structure between 1.8 GeV/c^ and 3.0 GeV/c^ with a 
peak cross section of 10-30 nb near 2.2 GeV/c^, while 
the one-pion-exchange model predicted an enhance- 
ment near threshold around 1.6 GeV/c^ with a peak cross 
section of 13 nb using a preferred value of the slope pa- 
rameter. Both the pea k position and the peak height pre- 
dicted in [2^ and 24 1 disagree with our measurements. 
The products of the two-photon decay width and branch- 
ing fraction of the rjc, XcO and Xc2 to cj0, 00 and ujtu are 
also measured assuming no interference. 
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